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The canonical view of the interactions between viruses and their microbial hosts presumes that 
changes in host and virus fate require the initiation of infection of a host by a virus. That is, 
first virus particles diffuse randomly outside of host cells, then the virus genome enters the target 
host cell, and only then do intracellular dynamics and regulation of virus and host cell fate unfold. 
Intracellular dynamics may lead to the death of the host cell and release of viruses, to the elimination 
of the virus genome through cellular defense mechanisms, or the integration of the virus genome 
with the host as a chromosomal or extra-chromosomal element. Here we revisit this canonical 
view, inspired by recent experimental findings of Bautista and colleagues (mBio, 2015) in which the 
majority of target host cells can be induced into a dormant state when exposed to either active or 
de-activated viruses, even when viruses are present at low relative titer. We propose that both the 
qualitative phenomena and the quantitative time-scales of dormancy induction can be reconciled 
given the hypothesis that cellular physiology can be altered by contact on the surface of host cells 
rather than strictly by infection. We develop a biophysical model of contact-mediated dynamics 
involving virus particles and target cells. We show how in this model virus particles can catalyze 
- extracellularly - cellular transformations amongst many cells, even if they ultimately infect only 
one (or none). We discuss implications of the present biophysical model relevant to the study of 
virus-microbe interactions more generally. 


I. INTRODUCTION 

Dormancy is ubiquitous in microbial systems. One 
prominent example is that of bacterial persistence, in 
which cells undergo a rapid physiological change char¬ 
acterized by slowed or even halted growth as well as 
decreased susceptibility to antibiotics w- Another 
example is that of microbial “seed banks” [4j|, in which 
individual microbes undergo a long period of stasis in 
which growth is halted. These cells can then restart 
growth given changes in environmental conditions. A 
third example is the starvation-dependent division of 
division of the bacteria Sinorhizobium meliloti into dis¬ 
tinct daughter phenotypes, one that is better suited to 
long-term starvation [5|. Long-standing theory |6] and 
recent experiments [7] support the consensus that vary¬ 
ing phenotypes - whether by stochastic bet hedging or 
via phenotypic plasticity - is an evolutionarily favorable 
strategy in the face of uncertainty in environmental selec¬ 
tion pressures 0. One of the selection pressure that 
microbial cells face is the possibility of infection and lysis 
by viruses. Indeed, it has been suggested that E. coli may 
enter the persistence state as a route to diminish, tem¬ 
porarily, the ability of viruses to eliminate a local pop¬ 
ulation [5]. Nonetheless, the transition between active 
growth and persistence was not proposed to be a func¬ 
tion of virus-host interactions. 

Here, we propose a biophysical model of virus- 
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induced dormancy of microbial host cells. Our model 
is inspired by recent empirical findings by Bautista and 
colleagues m of the interactions between the archaeon 
Sulfolobus islandicus and the dsDNA fusellovirus Sulfu- 
lobus spindle shaped virus (SSV9). S. islandicus is a 
globally distributed archaeon, commonly found in hot 
spring ecosystems. S. islandicus is also a model sys¬ 
tem for studying the eco-evolutionary basis for diver¬ 
sity in archaea [TTHl5j . A recent study of the inter¬ 
actions between the host strain S.islandicus RJW002 
and the virus SSV9 were recently shown to lead to a 
novel population-level outcome: “challenge of RJW002 
with SSV9 induced a population-wide stasis or dorman¬ 
cy response, where the majority of cells are viable but 
not actively growing” m- Dormant cells appear “emp¬ 
ty” without coherent intracellular structure in contrast 
to normal cells. Dormant cells can then reorganize and 
revert to actively growing cells. In the experiment, virus¬ 
es were introduced at low concentrations relative to that 
of hosts. Yet, after 24 hrs, nearly 100% of cells were 
classified as dormant m- In other words: there was a 
large-scale amplification in the number of dormant cells 
at the end of the 24 hr period vs. the number of viruses 
at the start of the experiment. Further, in a follow-up 
experiment, nearly 100% of cells initiated dormancy even 
when the host was exposed to de-activated viruses at low 
relative concentration. 

Bautista and colleagues [lOj highlight the potential role 
of dormancy as a strategy to survive viral infection and 
lysis. These experiments also raise the possibility that 
contact between viruses particles and host surface may be 
sufficient to initiate a large-scale physiological response, 
both at the cellular and population-scales. The model we 
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FIG. 1: Schematic representation of a biophysical model of virus-host interactions including susceptible microbial cells (S), 
complexes of a cell with a virus particle (C), infected cells (/), dormant cells ( D ) and free viruses (V). All populations are 
tracked in terms of densities of particles/ml. 


propose here focuses on the early dynamics of virus-host 
interactions. We assume that viruses can “contact” host 
cells reversibly. If most contacts between virus particles 
and host cells do not lead to infection, then there can be 
a broad dynamic regime in which nearly all of the target 
host cells enter dormancy even when there are far few¬ 
er viruses than hosts. We reach this conclusion by first 
proposing a biophysical model of contact-mediated dor¬ 
mancy. Next, we solve the model and identify critical 
transitions between distinct qualitative regimes, includ¬ 
ing a broad regime of dormancy-enhancement. Final¬ 
ly, we apply the model to the experimental conditions 
underlying the recent study of RJW002 and SSV9 and 
show that a contact-model of dormancy is compatible 
with observations. We discuss the implications of the 
current work for updating the dominant virus-host infec¬ 
tion paradigm to include renewed focus on the role of 
virus particles in stimulating changes amongst microbes 
and microbial populations. 


II. RESULTS 

A. Biophysical model of contact-mediated 
dormancy of hosts by viruses 


We propose a nonlinear dynamics model of virus-host 
interactions (see Figure [l]). This model describes the ear¬ 
ly dynamics of interactions involving viruses and hosts 


and the initiation of either dormancy or an active infec¬ 
tion. Consider an environment containing susceptible 
cells, S, and free virus particles, V. Free viruses can 
contact cells forming a complex, C, given a diffusion- 
limited contact rate of k + cells/(mhhrs). The use of the 
term “complex” suggests an analogy to models of enzyme 
kinetics. The complex is reversible. The disassociation 
rate is fc_ cells/(ml-hrs) and the infection rate is kf hrs -1 . 
If disassociation takes place, then the virus is released 
back into the environment. We assume that disassoci¬ 
ation may also induce a cellular transformation leading 
to dormancy with probability p. If infection takes place, 
then the virus genome enters the host cell, leading to an 
actively infected cell. The dynamics of this model can be 
written as: 
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The system can be reduced in complexity. First there 
is a constraint that So = S(t)+C(t)+D(t)+I(t), because 
the dynamics in the model track the transformation of an 
initial population of So susceptible cells into four differ¬ 
ent states: susceptible, complex, dormant, and infected. 
There is another constraint that Vq = C(t)+I(t) + V(t), 
because the dynamics track the transformation of an ini¬ 
tial population of Vq virus genomes into three different 
states: in free viruses, temporarily bound with cells, and 
injected into hosts. Finally, when contact occurs rapidly, 
then we can use a standard assumption in enzyme kinet¬ 
ics theory and presume that the concentration of C rapid¬ 
ly equilibriates (for example, see the Appendix of jT5]). 
This is a standard approach to analyzing models char¬ 
acterized by fast-slow dynamics. Here, we assume that 
the change in C is relatively fast when compared to oth¬ 
er state variables. In the fast limit, then C(t) = k + ^ f ■ 
This approximation is referred to as quasi-steady state 
approximation (QSSA). Substituting the QSSA equilib¬ 
rium for the concentration of the complex yields the fol¬ 
lowing reduced system: 


d S 
dt 
dV_ 
dt 


= -SV 
= -SV 


kfk + + pk_k + \ 
k- + k f ) 
k +kf \ 
k- + kf ) 


( 2 ) 


We can then identify the following control parameters: 
the conditional probability of infection given contact, q = 
k k + kf , the effective adsorption rate, <j> = qk+, and the 
ratio of dormancy induction to infection, 5 = . Using 

these control parameters, we can rewrite the model as: 


_ = —<t>SV (1 + 5) (3) 

dV 

- = -*SV (4) 


This model can be interpreted as follows. The density 
of susceptible hosts decreases at a rate proportional to 
the densities of virus and susceptible host populations 
(Eq. ©). The proportionality constant is </>, the adsorp¬ 
tion rate, multiplied by an enhancement factor of (1 + 5), 
where 5 is the number of dormant cells induced for each 
infected cell produced. The enhancement factor arises 
due to the fact that susceptible hosts can become infected 
or enter dormancy due to interactions with viruses. The 
number of free viruses decreases at a rate proportional 
to the densities of virus and susceptible host populations 
(Eq. 0). The proportionality constant in that case is </>, 
the adsorption rate, because that is the means by which 
free viruses are removed from the medium. 

The model does not include birth and death of hosts, 
the lysis of hosts by viruses, nor the decay of virus par¬ 
ticles. As such, the model describes early dynamics 
of virus-host interactions. This focus is in contrast to 
models of virus and host dynamics mediated by density- 
dependent infection and lysis EM]. 
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TABLE I: Asymptotic densities of state variables given the 
control parameter fi = So — (1 + 5)Vo. The conditions U < 0 
and > 0 represents the host-depletion and virus-depletion 
limits. See the text for more details. 


B. Qualitative regimes of dormancy induction 


Eqs. ([3|-0 can be solved analytically (see 
Appendix]A|, yielding: 


S(t) 

V(t) 


n 

1-(1 + 5)Moe-^ t 
M 0 fle~^ nt 
1 - (1 + 5)Moe~^ nt 


(5) 

( 6 ) 


where Q = So — (1 + S)V 0 and M. 0 = Vq/Sq■ These 
solutions hold so long as O / 0. When fl = 0 then 


S(t) = 

So 

(7) 

1 + Soffit 

V{t) = 

Vo 

(8) 

1 + Vo4>{ 1 + 5)i 


The system dynamics have qualitatively different behav¬ 
iors for Q > 0 and for Q < 0 (see Table [I]). Therefore 
Q acts as a critical parameter, both in a biological and 
dynamical systems sense. 

Recall that (1 + 5)+o is the maximum number of hosts 
that can be infected or enter dormancy as a result of 
interactions with viruses. Therefore, when So > (1 + 5) Vo 
then there are enough hosts for all viruses to infect cells 
(Vo in total) and to catalyze 5 hosts per infected cell to 
enter dormancy (5Vo in total). This is the case when 
Q = 5b — (1 + 5)Vo is positive. In this limit, all viruses 
infect a cell, while some hosts remain uninfected. The 
condition Q > 0 represents the “virus-depletion” limit. 
In contrast, when So < (1 + 5)Vo there are not enough 
hosts for all viruses to infect cells and to catalyze 5 hosts 
per infected cell to enter dormancy. This is the case when 
VI = Sq — (1 + 5)Vo is negative. In this limit, all hosts 
are either infected or enter dormancy, while the viruses 
remain in the system. The condition Q < 0 represents the 
“host-depletion” limit. The condition Q = 0 represents 
the critical point dividing these two dynamical regimes 
(see Figure [2]). 
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FIG. 2: Dynamics of susceptible hosts, S(t ), and free viruses, 
V(t), for the reduced model in Eq. (3|- Q in three regimes: 
fi > 0, fi = 0 and Q < 0 (top, middle and bottom, respec¬ 
tively). Common parameters for the dynamics are <j> = ICC' 
cells/(ml-hrs), S = 4.5 and So = 10' cells/ml. The initial 
virus density is Vo = 1.8 x 10 5 , 1.8 x 10 6 and 3.5 x 10 6 virus¬ 
es/ml in the top, middle and bottom panels respectively. 


C. Viruses can induce nearly all hosts to enter 
dormancy, even when the virus-host ratio is far less 
than one 

Traditional analysis of virus-host interactions presup¬ 
poses that entrance, through injection or other means, of 
virus genomes into a host is required for virus-mediated 
modification of host cell physiology. Here, as in tradi¬ 
tional models, Vq represents the upper limit to the num¬ 
ber of hosts infected by viruses (see Table [i]). This limit 
holds when restricting attention to short-term dynamics 
before replication and lysis which releases more viruses 
that can initiate subsequent infections. However, in the 
present model, hosts can also undergo contact-mediated 
dormancy. When the ratio of viruses to hosts is small, 
i.e., the MOI is Mo = Vo /So -C 1, we find an unexpect¬ 
ed outcome: nearly all of the hosts can enter dormancy 
even when there are far fewer viruses than hosts. In the 
host-depletion regime, H < 0, then D a0 = If 5 3> 1 
then Doo —> So- This condition holds so long as the rela¬ 
tive rates of unbinding and dormancy are high relative to 
infection and there are enough viruses. The critical virus 


density depends on 6 and is equal to So/(l+<5). For virus 
densities above this value, then the dormant cell fraction 
will have reached its maximum because the system moves 
from being host-depleted to virus-depleted. The compar¬ 
ison of the asymptotic dormant cell fraction and infected 
cell fraction are shown in Figure [3] As is apparent, more 
cells become dormant and infected with increasing titer. 
Yet, the balance of dormancy or infected cell fates shift 
with increases in S. As 6 increases, then many dormant 
cells are initiated for each infected cell, whereas when <5 
decreases, then very few dormant cells are initiated for 
each infected cell. 


D. Dynamics of dormancy induction in the 
archaeon Sulfolobus islandicus 

We apply our biophysical model of contact-mediated 
dormancy to a recent empirical study of interactions 
between the archaeon S. islandicus and the dsDNA 
fusellovirus Sulfulobus spindle shaped virus (SSV9). In 
the experiment, viruses were introduced at relatively low 
concentrations compared to that of hosts. The ratio of 
viruses to hosts, M was estimated to range between 0.01, 
given plaque-forming unit counts, and 0.01, given quan¬ 
titative PCR counts. In this experiment nearly 100% 
of host cells entered dormancy. Hence, there was a 10- 
fold to 100-fold increase in the conversion of host cells 
into a dormant state. SSV9 was then exposed to UV 
light to de-activate the virus population. The subsequent 
conversion estimates into dormant cells were statistically 
unchanged. We interpret this result to mean that con¬ 
tact between viruses and host may be sufficient to induce 
a population-wide dormancy response even when viruses 
are present at densities far below that of target host cells. 
These qualitative results are the basis for our quantita¬ 
tive parameterization and analysis of the model. 

The governing parameters of the biophysical model are 
</>, the adsorption rate, and 6, the ratio of dormancy 
induction to infection. Bautista and colleagues estimat¬ 
ed <f>original to be 8.4xl0~ n ml/min based on the decay 
of plaque-forming units. They estimated the adsorption 

rate using the formula Original = where V 0 is 

the original titer of viruses, V t is the titer at time t and N 0 
is the original titer of hosts. Conventional estimates are 

Jog ^-2. 

that 4> = , hence we downward adjust the adsorp¬ 

tion rate to be <f> = 2.2 x ICC 9 ml/hr (note the change 
in units). The effective adsorption rate is a combination 
of the process of diffusion-limited contact and success¬ 
ful infection. We estimate the diffusion-limited contact 
rate based on physical principles standard in the study 
of virus-host interactions mm-- 

k- l«4x ICC 9 — ml/hr (9) 

Th 

where r v is the effective radius of the virus, r^ is the 
effective radius of the host, where the prefactor is appro- 
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FIG. 3: Asymptotic fraction of dormant and infected cells resulting from virus-host dynamics. (Left) Final dormant cell fraction 
as a function of the initial ratio of viruses to hosts, Vo/So- (Right) Final infected cell fraction as a function of the initial ratio of 
viruses to hosts, Vo/So- As shown, the final cell state fractions depend on 5 with qualitatively different responses as a function 
of the control parameter S2. The fraction of dormant cells and infected cells increases linearly with Vo/So so long as > 0. 
When Q < 0 then the fraction of dormant and infected cells is a constant, irrespective of Vo/So- Increases in 8 lead to a relative 
increase in the proportion of dormant vs. infected cells. Parameters are otherwise the same as in Figure [2] 


priate for interactions taking place at room temperature 
(293 degrees K) and in a medium with the viscosity of 
water. Assuming 77, = 1 /1 m and r v = 0.04 /im then 
we predict k + = 10“' ml/hr. The ratio of the diffusion- 
limited contact rate expected from first principles, 
and the realized adsorption rate measured in the exper¬ 
iment, 4> can be used to estimate q — <f>/k + = 0.022. In 
this limit, then q w kf/k-, such that 6 = p/q. 

We can explore the predicted fraction of dormant cells 
as a function of S. First, we analyze the ratio of dormant 
cells to initial viruses given variation in q and p (Fig¬ 
ured-left). For q <C 1 there is a broad range of values of 
p such that Doo/Vq 3> 1. This result means that many 
hosts cells could enter dormancy given exposure to low 
relative titer of viruses. Yet, the value of p, and there¬ 
fore of 8, remains a free parameter given the experimental 
tests conducted in the system. The maximum value of <5 
is when p -A 1, i.e., when the conditional probability of 
inducing dormancy given a reversible contact approach¬ 
es 1. In this case, <5 -A 1/0.022 = 45. Lower values of 8 
are possible when p —> 0, i.e., when the conditional prob¬ 
ability of inducing dormancy given a reversible contact 
approaches 0. In this limit, <5 —> 0. The experimen¬ 
tal finding of a range between 8 = 10 to 8 = 100 fold 
enhancement is consistent with our finding of 8 = 45 and 
a model in which p —► 1 (see Figured-right). 

We provide another evaluation of the model by con¬ 
sidering the timescale over which contact-mediated dor¬ 
mancy should take place. The appropriate time-scale is 
predicted to be r c = ■ The approximate cell den¬ 

sity for experiments in [TO] was Sq ~ 2.5 x 10 8 ml” 1 . We 


assume that viruses were present at Vo « O.OlS’o- Using 
these values we estimate r c ~ 4 hrs. Hence, we predict a 
characteristic time-scale for conversion of 64% of hosts, 
corresponding to a one-log drop in susceptible host den¬ 
sity, in a time period of 4 hrs and to conversion of 87% of 
hosts, corresponding to a two-log drop in susceptible host 
density, in a time period of 8 hrs. We view this time-scale 
analysis to be another confirmation of the model, given 
that even if the hosts initiate dormancy after contact, 
the dormancy “phenotype” is likely to be delayed given 
the re-organization of intracellular dynamics. In sum¬ 
mary, non-infectious contacts could happen sufficiently 
frequently so as to rapidly induce dormancy on relevant 
time-scales of experimental observations. 


III. DISCUSSION 

We have proposed a biophysical model of host- 
dormancy initiated by contact with viruses. The mod¬ 
el explicitly accounts for the possibility that viruses can 
contact host cells reversibly. Reversible contacts may, 
with some frequency, lead to induction of dormancy. 
Such contact-mediated dormancy at the cellular scale 
can be evident at the population-scale in certain lim¬ 
its. In particular, we predict a critical transition to a 
regime in which the vast majority of cells become dor¬ 
mant even if the initial ratio of viruses to hosts is quite 
small. This regime is found to be robust to a broad range 
of biophysically-relevant parameters. 

The inspiration for the model was a recent series of 
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FIG. 4: Ratio of the concentration of dormant cells at the end of the dynamics to the concentration of viruses at the start 
of the dynamics. The intensity values denote the ratio, D a0 /Vq. The solid line demarcates the boundary between ratios that 
exceed one (upper-left portion) and those that are less than one (remainder). (Left) Population dormancy enhancement as a 
function of the probability of infection given contact, q, and the probability of dormancy-initiation given reversal of contact, 
p. Here, So = 2.5 x 10 8 cells/ml and Vo = So/100. (Right) Population dormancy enhancement as a function of the initial 
virus-host ratio, Vo/So, and the cell fate ratio, 5 ~ p/q. Key point: the number of cells that enter dormancy per virus can be 
much greater than one, even if the initial virus-host ratio is much smaller than 1. 


findings that the majority of an archaeal population 
could enter a stasis-like “dormancy” in less than 24 hrs 
after exposure to a relatively small number of viruses m- 
The same effect was observed whether active or inacti¬ 
vated viruses were utilized. This experimental finding 
suggests the possibility that contact between virus par¬ 
ticles and host surfaces induce a transformation in host 
phenotype. Dormant cells were unlikely to be infected 
and lysed by viruses. However, such dormancy comes at 
a cost, as residing in a dormant state for too long can lead 
to loss of cell viability and cell death pm- This transfor¬ 
mation may represent a form of phenotypic plasticity on 
the part of hosts. Further work would be required to 
evaluate whether dormancy could be initiated indepen¬ 
dent of virus-contact, which would represent a form of 
bet hedging. 

We fit our model to the experimental host-virus sys¬ 
tem, leaving only one free parameter: the conditional 
probability of dormancy initiation upon a reversible con¬ 
tact. We predict that whenever this conditional probabil¬ 
ity is sufficiently high, then large-scale initiation of dor¬ 
mancy can occur even when very few cells are infected. 
Based on our fits, we predict rapid initiation of dorman¬ 
cy can take place on a time-scale of 4-8 hrs, sufficiently 
fast so as to identify a dormancy phenotype amongst the 
majority of the host population in the 12-24 hrs period 
as observed. Moreover, the same model fits predicts the 
potential for a 50-fold enhancement in dormant cells with 
respect to viruses. Experiments observe higher enhance¬ 
ment ratios ranging from 10- to 100-fold m- The uncer¬ 
tainty is due, in part, to challenges in quantifying virus 


titer. Yet, there are other challenges, including poten¬ 
tial errors in the estimation of the adsorption rate and 
approximation of the maximum contact rate. There may 
also be additional mechanisms of relevance, e.g., cell¬ 
cell communication, as a means to amplify a small viral 
contact “signal” or additional spatial structure of the 
cells and viruses in the environment not accounted for 
in the present model. Nonetheless, the current biophys¬ 
ical model provides an explanation for the qualitative 
features of dormancy-enhancement and the time-scale of 
the effect. We suggests that the present baseline model 
serve as the basis for future detailed investigations of the 
phenomena moving forward. 

A number of issues remain to link the proposed early- 
time dynamics with long-term dynamics. First, nonlin¬ 
ear feedbacks are likely to arise in this system due to 
the infection and release of viruses. These viruses are 
themselves metastable, and so incorporating the decay 
of infectious virus particles will also need to be consid¬ 
ered. Second, here we assume that viruses cannot infect 
dormant cells. The interaction between viruses and dor¬ 
mant cells is not fully elucidated. Finally, it is known that 
intracellular interactions of S. islandicus and its viruses 
are mediated, in part, by the CRISPR/Cas immune sys¬ 
tem. The CRISPR/Cas immune system is ubiquitous 
in bacteria and archaea. CRISPR/Cas enable host cells 
to target and degrade foreign genetic elements, including 
viruses [2TlT - f2o| . CRISPR-mediated interactions can lead, 
over time, to the diversification of the host as it obtains 
new immune elements from the virus and to the diversi¬ 
fication of the virus [2BH251 . Linking early- to long-term 
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dynamics will also need to confront the potential diversi¬ 
fication of communities arising due to contact-mediated 
and infection-mediated dynamics. 

In summary, dormancy is a feature of organisms span¬ 
ning animals to plants to microbes. The evolution of dor¬ 
mancy has long been thought to represent a way to max¬ 
imize long-term fitness in an uncertain environment [B]. 
For microbes, part of the uncertainty in their fitness 
stems from the possibility that they may be infected and 
lysed by a virus. Here, we find that a biophysical mecha¬ 
nism of context-initiated dormancy can explain observa¬ 
tions of rapid and large-scale conversion of a host archaeal 
population into a dormant state by a relatively small 
number of viruses. We suggest that detailed investiga¬ 
tions of contact of hosts by viruses is likely to yield new 
biological surprises. In turn, new theoretical approaches 
are needed to consider the integration of “fast” dynamics 
at contact-scales with the long-term nonlinear feedbacks 
arising from the effects of physiological transformations 
and infection on host and virus populations. 
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Appendix A: Analytical solution of the 
contact-mediated dynamics model 

The following reduced model introduced in Eqs §-0 repre¬ 
sents dynamics of susceptible hosts, S(t ), and free viruses, V(t): 

-</>SV(l + 6) (Al) 

-(/>SV (A2) 


d5 
d t 
dV 
d t 


This model can be solved analytically by observing that: 


**-(! + *) 
dV y J 

(A3) 

such that 


s(t) = (i + <5)vx*) + n 

(A4) 

where the integration constant, can be solved 

conditions that 5(0) = So and V(0) = Vo: 

using the initial 

O = 5 0 -(l + <5)Vb 

(A5) 


This two-dimensional system can be reduce to a one-dimensional 
system by substituting S(t ) as in Eq. (A4) , yielding: 

=-^y((i + 5)y + n) (A6) 

When £7 0, then Eq. ( | \S[ i can be solved by separation of 

variables yielding: 


V(t) 


Motte-'f’™ 

1 - (1 + 5)M 0 e~^ m 


(A7) 


where Ado is the initial population-level multiplicity of infection, 
i.e., the ratio Ado — Vo/S'o- The solution for V(t) is the basis for a 
complete description of the dynamics under the Quasi-steady state 
approximation: 


S(t) 

v(t) 


n 

1 — (1 + S)Moe~ c l ,nt 
_M 0 ne-^ nt 
1 - (1 + (5)Aloe“^ nt 


(AS) 


where, in addition: 


o( t) = * S(t)V(t) 

k f 

m = Vo - C(i) - V(t) 

D(t) = (So - S(t)) - (Vo - V(t)) 


(A9) 


When £7 = 0, then Eq. (| A6|) can be solved by integration: 


S(t) 

V(t) 


So 

1 + Soffit 
Vo 

1 + Voifi( 1 + <5 )t 


(A10) 


recalling that when £7 = 0 then (1 + <5)Vo — So. 
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